Smooth Hybrid Inflation 

G. Lazarides and C. Panagiotakopoulos 

Physics Division 
School of Technology 
University of Thessaloniki 
Thessaloniki, Greece 



Abstract 

We propose a variant of hybrid inflation which is applicable in a wide 
class of supersymmetric grand unified models. The observed temperature 
perturbations of cosmic background radiation can be reproduced with natural 
values of the parameters of the theory and a grand unification scale which 
is consistent with the unification of the minimal supersymmetric standard 
model gauge couplings as measured at LEP. The termination of inflation is 
smooth and does not produce any topological defects. Finally, we present a 
specific supersymmetric model where our smooth hybrid inflationary scenario 
is realized. 
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The usual realizations of the new and the chaotic inflationary scenarios*- 1 - 1 
invoke a very weakly coupled gauge singlet scalar field known as inflaton. The 
extremely weak couplings of this field, which are necessitated by the smallness 
of the observed temperature fluctuations in the cosmic background radiation 
(CBR) are certainly unsatisfactory because they create an extra 'hierarchy' 
or 'naturalness' problem. Recently, Linde*^ has proposed a new realization 
of the chaotic inflationary scenario based on a coupled system of two scalar 
fields one of which may not be a gauge singlet. The great advantage of 
this hybrid inflationary model is that it produces the observed temperature 
fluctuations in the CBR with natural values of the coupling constants. How- 
ever, inflation terminates abruptly and is followed by a 'waterfall' regime 
during which topological defects can be easily produced. Consequently, in 
the cases where these defects include superheavy magnetic monopoles or do- 
main walls we end up with an unacceptable cosmology. Of course, this is not 
the case with the original hybrid inflationary model*- 2 - 1 . Hybrid inflation has 
been adjusted*- 3 -* so that it becomes applicable in a wide class of supersym- 
metric (SUSY) grand unified theories (GUTs) based on semi-simple gauge 
groups. The magnitude of the temperature fluctuations in CBR turns out 
to be directly related to the grand unification scale M X - However, values 
of Mx consistent with the unification of the gauge couplings* 4 -' of the min- 
imal supersymmetric standard model (MSSM), which is favored by recent 
LEP data, tend to give values of the temperature fluctuations in CBR con- 
siderably higher than the values allowed by observations. Also, the scheme 
heavily relies on radiative corrections. 

In this paper we propose a variant of Linde's potential which can also be 
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derived in a wide class of SUSY GUTs based on semi-simple gauge groups. 
Instead of the renormalizable trilinear coupling in the superpotential we uti- 
lize the first non-renormalizable contribution. All our analysis is done at the 
tree-level and no radiative corrections are needed. Although we get only a 
slight variation of Linde's potential, the cosmological scenario obtained is 
drastically different. Already since the beginning of inflation, the system 
follows a particular valley of minima which leads to a particular point of the 
vacuum manifold. So our inflationary scenario does not lead to production 
of topological defects. Also the termination of inflation is not as abrupt as 
in the hybrid case. It is quite smooth and resembles more the cases of new 
or chaotic inflation. The main advantage of our smooth inflationary scenario 
is that the measured value of the temperature fluctuations of CBR can be 
reproduced with natural values of the parameters and with a GUT scale, 
Mx, consistent with the unification of the MSSM gauge couplings. It is also 
remarkable that the scale controlling the non-renormalizable terms in the su- 
perpotential turns out to be of order 10 18 GeV. The spectral index of density 
fluctuations is close to unity. 

We consider a SUSY GUT based on a (semi-simple) gauge group G of 
rank >5. G breaks spontaneously directly to the standard model (SM) group 
Gs at a scale Mx ~ 10 16 GeV. We assume that below Mx the only SM 
non-singlet states of the theory are the usual MSSM states in order for the 
successful MSSM predictions for sin 2 9 w and a s to be retained. The theory 
could also possess some global symmetries. The symmetry breaking of G to 
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Gs is obtained through a superpotential which includes the terms 

Here 0, is a conjugate pair of left-handed SM singlet superfields which 
belong to non-trivial representations of the gauge group G and reduce its 
rank by their vacuum expectation values (vevs), s is a gauge singlet left- 
handed superfield, \x is a superheavy mass scale related to Mx, whereas 
M is a mass scale of the order of the 'compactification' scale M c ~ 10 18 
GeV which controls the non-renormalizable terms in the superpotential of 
the theory. The superpotential terms in eq.(l) are the dominant couplings 
involving the superfields s, 0, consistent with a continuous R-symmetry 
under which W — > e W, s — > e l9 s, 00 — > 00 and a discrete symmetry under 
which 00 changes sign. Moreover, we assume that the presence of other SM 
singlets in the theory does not affect the superpotential in eq. (1) which is 
responsible for the non-zero vevs of the fields 0, 0. This, if not automatic, is 
very easily achieved in the case of semi-simple gauge groups. The potential 
obtained from W in eq.(l), in the supersymmetric limit, is 

^l/-ff + 4| S |« (W + \W) + D - terms, (2) 

where the scalar components of the superfields are denoted by the same sym- 
bols as the corresponding superfields. Vanishing of the D-terms is achieved 
with |0| = |0| (D-flatness condition). The supersymmetric vacuum 

< s >= 0, < >< >= ±(J,M, | < > | = | < > | (3) 

lies on the D-flat direction 0* = ±0. Restricting ourselves to this particular 
direction and performing appropriate gauge, discrete and R-transformations 



we can bring the complex s, 0, fields on the real axis, i.e, s = = = 
\x, where a and x are re& l scalar fields. The potential in eq.(2) then takes 
the form 



%-) = (^ 2 -i) 2 + Sl (4) 



c 4 )2 | xV 2 

16M 2> 16M 4 
and the supersymmetric minima correspond to 



| < X > I = 2(/iM) 1 / 2 , < a >= 0. (5) 

The mass acquired by the gauge bosons is M x = g (/iM) 1 / 2 , where g is the 
GUT gauge coupling. The vacuum manifold of the theory, which is obtained 
from the minimum in eq.(5) by performing all possible gauge and global 
transformations, may have non-trivial homotopical properties in which case 
the theory predicts the existence of topological defects. 

Let us now take a closer look at the potential in eq. (4). For any fixed 
value of a - , this potential, as a function of x 2 , has a local maximum at x 2 — 
and an absolute minimum lying at 



4/x 2 M 2 
3^ 



X 2 ^ o^^, f° r a » ( 6 ) 



and x 2 — 4/iM, for a 2 « fiM. The value of the potential along the maxima 
at x 2 — is constant, V max (x 2 = 0) = /i 4 , whereas along the valley of minima, 
for a 2 » fiM, is V min (a) ~ /i 4 [l - (2/27)(/i 2 M 2 /cr 4 )]. 

We are now ready to turn to the discussion of the cosmological evolution 
of this system. We assume that after 'compactification' at a cosmic time 
t c ~ (|^) 1/2 §t, where M P = 1.2 x 10 19 GeV is the Planck mass, the universe 
emerges with energy density of order M 4 . We then, as usual, follow the 
development of a region of size £ c ~ t c where the scalar fields x an d o happen 



to be almost uniform with |er| >> \x\- Under these circumstances, one can 
easily see that initially the last term in eq.(4) is the dominant contribution 
to the potential energy of the system. Thus, the initial equations of motion 
for the x an d ° fields read 

3(jV 



and 



a+3H&+ m^°> (8) 



with H being the Hubble parameter and dots denoting derivatives with re- 
spect to cosmic time. Assuming for the moment that a remains approx- 
imately constant, we see that the frequency of oscillation of the x field, 
which is of order Xm\ a \/M 2 (xm is the amplitude of the oscillations of x), 
is much greater that the Hubble parameter H ~ (7r/6) 1 / 2 (x^ t |o"|/MpM 2 ), 
for Xm « Mp. So x initially performs damped oscillations over the maxi- 
mum at x — 0- For |er| >> Mp, the 'frequency' of oscillation of the a field, 
which is of order x m /M 2 , is much smaller than H . Thus, our hypothesis 
that a does not change much is justified. When the amplitude of the x field 
drops to about (/x 2 iW 2 /|cr|) 1,/3 , the /i 4 term dominates the potential in eq. 
(4) and the Hubble parameter becomes approximately constant and equal 
to H = (87v/3) 1/2 fi 2 /M P and remains so thereafter till the end of inflation. 
However, the form of the eqs.(7) and (8) still holds till Xm drops to about 
fxM/\a\. The 'frequency' of the a field remains much smaller than H, for 
(/i 2 M 2 /M) 1/3 £ Xm Z tiM/\(r\ (H » M P ) and eq.(8) reduces to 
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As is easily shown, a again remains almost constant within one expansion 
time, i.e., \Aa/a\ « x<L m p/ ' < OWN) 2 « l - Th e frequency of 
the x field ~ XmW\/M 2 remains greater than H, for Xm £ [iM / (M P \a\) 1 / 2 > 
fxM/\a\ and x still performs damped oscillations. For smaller values of the x 
field (jj,M/(Mp\a\Y^ 2 Z x ~ A 4 ^/! "!)? the Hubble parameter overtakes and 
X enters into a 'slow roll-over' regime controlled by the equation 

3/f * + W" ' (10) 

Eqs.(9) and (10) then imply that 3a 2 — x 2 — constant. Thus, x keeps drop- 
ping without any essential change in a and, as one can deduce from eq.(10), 
reaches a value ~ \iMj\a\ in a cosmic time interval At ~ (a / Mp) 2 H~ l » 
if -1 . After that, the x 4 term in eq.(4) comes into play and the equation of 
motion of the x field becomes 

5x + 3H5x + m 2 x (a)5 X ^0, (11) 

where 5x is the deviation of x from the minima in eq.(6) and m^a) ~ 
4///3<t 2 is the a -dependent mass squared of the x field along the valley of 
these minima. For |<r| >> M P , m x (a) « H and eq. (11) reduces to 

3H5x + m 2 x (a)5x ^ 0. (12) 
Thus, x reaches the minimum in a time interval 

M ~ 67r( ^ )2//_1 (13) 

within which the a field stays essentially unaltered as one can see from eq.(9). 
To summarize, so far, we have seen that within a time interval given in eq.(13) 



the x field falls into the valley of minima in eq. (6) and relaxes at the bottom 
of this valley whereas the a field still remains unchanged and much greater 
than M P . Of course, this does not automatically mean that x follows the 
valley of minima at subsequent times but, as we shall soon see, it turns out 
that it actually does. 

One could also reach the same conclusion starting with field values con- 
siderably smaller than Mp. This requires that the initial energy density be 
~ /i 4 for a to be slowly varying. Consider the situation where the initial 
value of a is close to M. It is then easy to show that, for initial values of x 
much smaller than (4\Z37r ^9 1 ) 1 / 3 ( < uM) 1 / 2 but still greater than the value 



M P 

given in eq.(6), the x field performs damped oscillations and ends up at the 
bottom of the valley of minima with a remaining essentially unaltered. 

To discuss the further evolution of the system let us, for the moment, 
suppose that Xi a ^ subsequent times, follows the valley of minima in eq. (6). 
Eq. (9) then reads 

3^ + 5^0, (M) 



which gives 

,2/iMAfp j 



a 



' /3JV ' /V) ' <15) 
where N(a) = HAt is the number of e-foldings from the moment at which 

the a field has the value a till the end of inflation. Inflation ends when the 

'frequency' of the a field becomes smaller than 3H/2, i.e., at 

^--(v^' 1 ' 3 ^) 1 ' 2 ^^ 1 ' 2 (16) 

For \a\ > M P , the evolution of x is governed by eq.(12) and the time needed 
to reach its minimun (see eq.(13)) is much smaller than A^(o")iJ _1 from 
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eq.(15) which is the time required for o to change significantly. Also the quan- 
tum fluctuations of x in de Sitter space, 5x — H/2n, are not strong enough to 
surpass the barrier at x = if cr << \^2tt(M/ fi)M P . For a < M P , m x > H, 
the quantum fluctuations of x stop and we go back to eq.(ll) which de- 
scribes damped oscillations about the minimum in eq.(6) with damping time 
H^ 1 much smaller than N(a)H~ l . Thus, we conclude that our starting as- 
sumption that x follows the valley of minima is correct. Note that, due 
to supersymmetry breaking, the x an d cr fields acquire additional masses of 
order M s ~ 1 TeV. It is easy to check that this mass contribution is never im- 
portant for the evolution of the x field. However, it plays a role for the a field 
if its initial value happens to be greater than a s = (2/y / 6)(M/M s ) 1 / 3 / u. This 
introduces only unimportant complications to the above discussion which we 
will skip. In any case, as soon as a drops below a s its SUSY breaking mass 
becomes subdominant. 

The contribution of the scalar metric perturbation to the microwave back- 
ground quadrupole anisotropy (scalar Sachs- Wolfe effect) is given^ 5 ) by 
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w 32tTi /2 V 3 / 2 7T 1/2 a , 

{ T )S ~ { 45 ) M P (dV/day k ~ H l 10 J M P M^ H) K } 

where the right-hand side is evaluated at the value of the a field where 
the lenght scale fc -1 , which corresponds to the present horizon size, crossed 
outside the de Sitter horizon during inflation. The derivative (dV/da) is, of 
course, calculated on the valley of minima of eq.(6) and substituting a in 
terms of N(a) from eq.(15) and fx in terms of M X) we obtain 

(^) 5 ~ 5-^ 2 6 1 / 3 n' 1 / 3 N 5 H /6 g~ 1 °/ 3 M 1 x 0/3 M p 4/3 M- 2 , (18) 
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where Nh is the number of e-foldings of the present horizon size during 
inflation. Taking N H = 60, M x = 2 x 10 16 GeV, g = 0.7 (consistent with 
the MSSM unification) and (AT/T) - 5x 10~ 6 from COBE, we get M ~ 
9.4 x 10 17 GeV and /i ~ 8.7 x 10 14 GeV. We ignored the gravitational wave 
contribution (^t-)t to the quadrupole anisotropy since it turns out to be 
utterly negligible relative to the scalar component: 

AT V 1 / 2 

(— ) r ~ 0.78^ -4.1 xKT 9 . (19) 

„-i 



The amplitude of the density fluctuations on a given length scale k as this 



scale crosses inside the post inflationary horizon is proportional to N^-i 



Ntf 6 [k 1 (Mpc) 1 10 4 ] 5 / 6Ar » , where N k -i is the number of e-foldings of this scale 
during inflation (see eq.(18)). The spectral index is then given by 

n = 1 ~ 0.97 (for N H ~ 60), (20) 

3N H 

which is very close to the Harrison- Zeldovich value (n = 1) and lies in the 
central range of values preferred by observations. We can also estimate the 
value of the a field at the end of inflation (see eq.(16)) a Q ~ 1.1 x 10 17 
GeV ~ 5AM X and its value when the present horizon size crossed outside 
the inflationary horizon (see eq.(15)) a H — (9N H /2) 1/6 a ~ 2.54<t ~ 2.7 x 
10 17 GeV. The total number of e-foldings during inflation (see eq.(15)) turns 
out to be greater than 10 12 , for initial \a\ > Mp, and SUSY breaking becomes 
irrelevant if the initial value of the a field does not exceed a s ~ 5.8M P (for 
M s ~ 1 TeV). 

The termination of inflation is not as abrupt as in the hybrid inflationary 
scenario where the slow roll-over regime is followed by the 'waterfall' . In 
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our case, after the end of inflation the a and x fields enter smoothly into an 
oscillatory phase about the global supersymmetric minimum of the potential 
in eq. (5) as in the 'new' inflationary scenario. The frequency of the oscillating 
fields is m a = m x = 2 v / 2(/i/M) 1 / 2 /i ~ 7.6 x 10 13 GeV. These fields should 
eventually decay into lighter particles and 'reheat' the universe so that its 
baryon asymmetry can be subsequently produced. We will postpone until 
later the discussion of these important issues because it depends strongly on 
the details of the particle physics model one adopts. 

Our inflationary scenario, in contrast to the hybrid one, produces no 
topological defects in the universe. This is because the a and x fields during 
inflation follow everywhere the same valley of minima which leads to a specific 
point of the vacuum manifold of the theory. As we discussed earlier, the x 
field stays close to the bottom of the valley and is unable to fluctuate over 
the maximum at x — 0- Also, for \a\ Z Mp, its mass is greater than H and, 
consequently, suffers no quantum fluctuations in de Sitter space. In hybrid 
inflation, the choice of the vacuum is made at the 'waterfall' after the end of 
inflation and, thus, topological defects can be easily produced. 

It is well-known that, as soon as one replaces global by local supersym- 
metry, the potential of the theory changes drastically and inflation becomes, 
in general, impossible. This is, to a large extent, due to the generation of a 
mass for the inflaton which is larger than H . In our model, one can easily 
check that no such mass is generated provided we employ the canonical form 
of the Kahler potential. This fact allows us to hope that, although, our dis- 
cussion is not expected to remain unaltered, the modifications necessitated 
by the inclusion of supergravity will not destroy the whole picture. 
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To construct an example of a SUSY model with the usual successful 
MSSM predictions of sin 2 9 w and a s where our smooth hybrid inflationary 
scenario is naturally realized, we consider a model with gauge symmetry 
group the subgroup G = SU(3) C x SU(2) L x SU(2) R x U(1) B -l x U(l) T of 
Eq and a global symmetry group C x B\ x 5 2 x R commuting with G. C is a 
Z2 matter parity, B\ and -B2 are Z\ symmetries and it! is a U(l) R -symmetry 
under which the superpotential W — > e ld W . The coupling constants of the 
five factors in G are assumed to be equal at a scale M c ~ 10 18 GeV, prob- 
ably related to a more fundamental theory Let q, q c , £, £ c , h, g, g c , N denote 
chiral superfields transforming under G as (3, 2,1, +1/3, +1), (3, 1, 2,-1/3, +1), 
(1,2,1,-1,+1), (1,1,2,+1,+1), (1,2,2,0,-2), (3,1,1,-2/3,-2), (3,l,l,+2/3,-2) and 
(1,1, 1,0, +4) respectively. Also let s denote a gauge singlet. The superfield 
content of the model consists of three q i: three qf, three £ i: three £\, one h, 
one h, one g, one g, one g c , one g c , one one £°, one £ a , one £ c , four q(3+ m ), 
four g m , four q , ( : 3 +m ), four q c m , one iV , one N a , three JVj, three iVj and one s 
(i = 1, 2, 3 and m = 1, 4). Under C, all superfields remain invariant except 
q iy q°, £i, £\. Under the generator of Bi, all superfields remain invariant except 
N a which changes sign and £ c - l £ c - l g- l g c which are multiplied by i. Under the 
generator of B 2 , all superfields remain invariant except iV which changes sign 
and h,g,g c which are multiplied by i. Finally under the U(l) R -symmetry, 
all superfields have charge +1/2 except N ,g,g c which have charge +1/3, s 
which has charge +1 and h,£ c ,£ c ,N which have charge 0. With the above 
charge assignments, the U(l) R -symmetry has color anomaly and therefore 
is a Peccei-Quinn symmetry solving the strong CP problem. 

The gauge group U{1) T breaks down at a scale M x through the vevs 
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of N ,N . The D-flat direction |iV | = |jV | is completely F-flat as well, 
because of the global symmetries. The large Yukawa couplings of the terms 
gg c N Q , h 2 N and gg c N Q could lead to a radiative change of the sign of the 
SUSY breaking mass squared term for one linear combination of the N Q 
and N Q bosons at a scale ~ Mx thereby generating a large vev for this 
combination along the D-and F-flat direction. The gauge group SU(2)r x 
U(1)b-l breaks down to U(l)y also at a scale M x through the vevs of the 
neutral components v c Q and v c Q in t c Q and l c . The relevant superpotential terms 
are the ones in eq. (1) with the role of <f> and played by v c Q and v c Q respectively. 
The vev of N Q breaks B\ down to a Z 2 subgroup, while N breaks B 2 down 
to a Z 2 subgroup. Besides, the pair of < N Q >, < N Q > breaks the U{1) R 
-symmetry down to a Z 2 discrete symmetry. Starting with N Q and N Q at 
their vevs we end up with the previously studied cosmological scenario. 

We assume that the allowed tree-level mass terms for the gauge non- 
singlets are all ~ Mx ■ Below the scale Mx the only Gs -non-singlet states 
are the ones of the MSSM. Ordinary light quarks and leptons are contained 
in the qi,qf,£i,£° fields all having negative matter parity. The electroweak 
Higgs doublets h^\h^ are the SU(2) R -doublet partners in h with mass 
~< NlNlv c v c > /M c 6 . An immediate consequence of this simple structure is 
the relation tan(3 ~ m t /mj among the vevs of h^ l \ and the 3rd generation 
quark masses. The g,g c pair could, in principle, mediate proton decay but 
the global symmetries stabilize the proton almost completely. 

Due to the appropriately chosen spectrum, the one-loop renormalization 
group equations above M x predict identical running for the gauge couplings 
of all the factors in G. Combining this fact with their assumed equality at 
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M c , we conclude that equality of the gauge couplings of Gs at the scale Mx is 
a justified boundary condition. The successful MSSM predictions for sin 2 9 w 
and a s follow immediately. 

Returning to the cosmologically more interesting neutral sector we ob- 
serve that we have exactly three right-handed neutrinos contained in the 
three £fs which acquire Majorana masses M u c ~ M 2 n 2 /Ml from the non- 
renormalizable terms vlv]v c {v c ^ jMl = 1,2,3). Therefore the v\ Ma- 
jorana masses are at most ~ 10 12 GeV. The simple structure of neutrino 
Dirac mass terms leads to the simple relation tanj3 ~ m^/m T for at 
least the 3rd generation. Thus, we expect a v T mass m VT ~ (tan(3m T ) 2 / M v c 
which makes the z/ T 's a significant component of the hot dark matter of the 
universe. The unbroken Z 2 matter parity C guarantees that the lowest su- 
persymmetric particle remains stable and contributes to the cold component 
of the dark matter together with the axion. 

The above mentioned terms providing the Majorana neutrino masses rep- 
resent the dominant couplings leading to inflaton decay with width V ~ 
j^M 3 fx 3 M~ 6 m, where m ~ 10 14 GeV is the inflaton mass. The corre- 
sponding 'reheat' temperature T r ~ \(TM P ) 1 / 2 ~ 10 11 GeV. The ufs pro- 
duced from the inflaton decay generate a lepton asymmetry^ 6 ) which later 
will be transformed into the observed baryon asymmetry through the non- 
perturbative baryon-and-lepton-number- violating effects in the standard model 
at temperatures ~ 1 TeV. 

The only topological defects predicted by this model are the domain walls 
associated with the spontaneous breaking of the discrete symmetries B\ and 
B 2 down to their Z 2 subgroups. These symmetries are, at the first place, 
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broken by the vevs of N a , N Q . Our assumption that these fields acquire 
their vevs before inflation then leads to the complete absence of catastrophic 
domain walls in the universe. 

In summary, in the context of SUSY GUTs based on semi-simple gauge 
groups, we constructed a smooth version of the hybrid inflationary scenario 
which can 'naturally' reproduce the observed temperature fluctuations in 
CBR with a GUT scale ~ 10 16 GeV consistent with the MSSM unification 
and a 'compactification' scale ~ 10 18 GeV. This inflationary scenario does not 
produce any topological defects. We also presented a specific model where 
our smooth hybrid inflation is realized. 
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